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Abstract. The evolution of the magnetic field is investigated for isolated 

\ as well as binary neutron stars. The overall nature of the field evolution is 

Q.^ ' seen to be similar for an initial crustal field and an expelled flux. The major 

ff^ i uncertainties of the present models of field evolution and the directions in 

' which further investigation are required are also discussed in detail. 

^ ' 
ON ; 

^ . 1. Introduction 

^: 

Q , There is no consensus regarding the generation of the magnetic field in neu- 

\—{ ' tron stars. The field could either be a fossil remnant from the progenitor 

c/3 . star in the form of Abrikosov fluxoids of the core proton super-conductor 

\ (Baym, Pethick &; Pines 1969, Ruderman 1972, Bhattacharya & Srinivasan 

^ ■ 1995). Or it could be generated after the formation of the neutron star in 

^ , which case the currents would be entirely confined to the solid crust (Bland- 

' ford, Applegate &: Hernquist 1983). Evidently, the nature of the evolution 

. I would depend very much on the internal field configuration. Observations 

and statistical analyses of existing pulsar data, nevertheless, indicate that 
significant decay of magnetic field is achieved only if the neutron star is a 
member of an interacting binary (Bailes 1989, Bhattacharya 1991, Hartman 
et al. 1997). 



The processes that are responsible for the field evolution in neutron stars in 
binaries are - a) expulsion of the magnetic flux from the super-conducting 
core during the phase of propeller spin-down, b) screening of the field by 
accreted matter and c) rapid ohmic decay of the crustal field in an accretion- 
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Figure 1. Evolution of the surface magnetic field in LMXBs. a) For initial crustal 

currents. With an wind accretion rate of M = 10^^'' Mq yr^^. Curves 1 to 5 correspond 
to initial current configurations centred at p = 10^'\ 10^^ ", 10^^, 10^^ ''', lO^^g cm~^. All 
curves correspond to an impurity content, Q = 0.0. A standard cooling has been assumed 
for the isolated phase. The wind and Roche-contact phases are plotted in expanded scales. 
The dashed and the solid curves correspond to accretion rates of M = 10^®, 10~^° M0 
yr~^in the Roche-contact phase, b) For an expelled flux. The dotted and the solid curves 
correspond to accretion rates of M = 10^^, 10^^" Mq yr^^in the Roche contact phase. 
The curves 1 to 5 correspond to Q = 0.0, 0.01, 0.02, 0.03 and 0.04 respectively. All curves 
correspond to a wind accretion rate of M = 1Q~^^ M© yr"'^. 



heated crust (for a review see Bhattacharya 2000). Diamagnetic screening 
of the field by accreted matter does not seem hkely to have any long-term 
effect (Konar 1997) and we shall exclude it from the present discussion. 
The other models invoke ohmic decay of the current loops for a perma- 
nent decrease in the field strength. In either case, the effect of accretion is 
two-fold. The heating reduces the electrical conductivity and consequently 
the ohmic decay time-scale inducing a faster decay. At the same time the 
material movement, caused by the deposition of matter on top of the crust, 
pushes the original current carrying layers into deeper and denser regions 
where the higher conductivity slows the decay down. The mass of the crust 
of a neutron star changes very little with a change in the total mass; ac- 
cretion therefore implies assimilation of the original crust into the super- 
conducting core. When the original current carrying regions undergo such 
assimilation, further decay is stopped altogether. Both the purely crustal 
model as well as the model assuming an expelled flux have been investi- 
gated by many authors (see Bhattacharya 2000). The important difference 
between our work and that of the other investigators lies in our assump- 
tion of a flux freezing upon the assimilation of the original current carrying 
layers into the super-conducting core. 
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TABLE I - Purely Crustal Field 



system 


final field and period 


comment 


isolated 
radio pulsars 


high field, long period 


no significant field decay 
in 10^ years 


HMXB 


high field, long period 


high-mass binary pulsars 
and solitary counterparts 


low field, long period 


not active as pulsars 


LMXB 


high field, long period 


high field low-mass 
binary pulsars 
and solitary counterparts 


low field, short period 


low field low-mass 
binary pulsars 
and solitary counterparts, 
millisecond pulsars 



2. Nature of Field Evolution 

2.1. GENERIC FEATURES 

The qualitative features of field evolution, as outlined below, are similar for 

a) an initial crustal field and b) an expelled flux. 

Pure Ohmic Decay in Isolated Neutron Stars (Konar 1997) 

1. A slow/fast cooling of neutron star implies a fast/slow decay; hence a 
low/high final field. 

2. Initial crustal currents concentrated at lower /higher densities gives rise 
to low/high final surface fields. 

3. Large impurity content makes the decay rapid and gives rise to smaller 
final fields. 

Accretion-Induced Field Decay in Accreting Neutron Stars 

(Konar k Bhattacharya 1997 - KBI) 

1. In an accreting neutron star the field undergoes an initial rapid decay, 
followed by slow down and an eventual freezing. 

2. A positive correlation between the rate of accretion and the final field 
strength is observed, giving rise to higher final saturation field strengths 
for higher rates of accretion. 

Magnetic Field and Spin Period - We have investigated the nature of 
the final 'magnetic field-spin period' combination. Our results agree well 
with the observations and are summarised in table-I (KBI, Konar & Bhat- 
tacharya 1999a - KBII). The nature of field evolution is similar for the 
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THE OBSERVATONAL PARADIGM 



NEUTRON STARS 



RADIO PULSARS 



X-RAY BINARIES 



LOW MASS 



HIGH MASS 



ISOLATED 



BINARY 



(MAGNET ARS?) 



LOW MASS 



HIGH MASS 



(MSP) 



Figure 2. The arrows indicate the expected evolutionary link between X-ray binaries 
and binary radio pulsars. 



model of spin-down induced flux expulsion (Konar & Bhattacharya 1999b 
- KBIII). Though there is one major difference as can be seen from fig.(l). 
To produce millisecond pulsars in LMXBs for an expelled flux large values 
of impurities, in the prior-to-accretion original crust, are required. But this 
would result in extremely small surface fields in old isolated pulsars. This 
is in complete contrast to a purely crustal model and expectation from 
statistical analyses of pulsar data. 

2.2. RANGES OF PHYSICAL PARAMETERS 

The paradigm of field evolution that have emerged out of observations, 
statistical analyses and theoretical expectations have been summarized in 
fig. (2) where the connection between the radio pulsars and their binary 
counterparts, namely the X-ray binaries is indicated. In table-II we indicate 
the constraints on various physical parameters in the field evolution models 
placed by the requirement to match observed properties in a variety of 
systems (KBII, KBIII). The parameters discussed here are - the density at 
which the initial crustal current distribution is located (pc), the impurity 
strength in the crust (Q), the duration of wind-accretion phase in different 
binary systems and the rate of accretion in the Roche-contact phase for 



LMXBs. 
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TABLE II - Constraints on Physical Parameters 



parameter 


model 


system 


requirement 


parameter range 


Po 


crustal 


HMXB 


high field 


high Po 


Q 


crustal 


isolated 

radio 

pulsar 


no decay over 
active pulsar 
life-time 


Q< 0.01 for 
standard cooling, 
Q< 0.05 for 
accelerated cooling 


expelled 
flux 


LMXB 


millisecond 
pulsar 
generation 


Q> 0.05 with 
wind accretion, 
Q >> 1 without 
wind accretion 


duration of 
wind accretion 


crustal 


HMXB 


high field 


short 


M in 

Roche-phase 


crustal 


LMXB 


high field 


Eddington rate 



3. Uncertainties and Future Directions 

The results and conclusions stated above suffer from a number of uncer- 
tainties regarding the micro-physics of the neutron star, as have been listed 
below. Moreover, a lot of the new theoretical results as well as the obser- 
vational facts have recently become available. In this section we mention 
some of the more important aspects that need to be incorporated in any 
future work on the evolution of the magnetic fields of neutron stars. 
Thermal Behaviour 

1. Isolated Phase - The present data can be made to fit scenarios with 
both a slow or an accelerated cooling. Therefore, it is not clear which is the 
correct cooling behaviour of an isolated neutron star. 

2. Accreting Phase - The crustal temperature corresponding to a given rate 
of accretion has not been determined with any degree of certainty. Also, 
the existing results are limited in their scope. 

3. Post-Accretion Phase - No calculation exists for the thermal behaviour 

of this phase at all. 

Transport Properties - Several factors affect the transport properties 
and hence both thermal and magnetic field evolution. Prominent among 
them is the change in the chemical composition due to a) accretion and 
b) spin-down. Recently it has been shown that the impurity content of the 
accreted crust for near-Eddington accretion rates could be extremely large 
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(Schatz et al. 1999). This, along with a temperature inversion near or be- 
yond the neutron drip (Brown 1999) might modify the transport properties 
significantly. Moreover the presence of dislocations, defects, non-spherical 
nuclei have so far not been taken into account in the calculation of trans- 
port properties. These are also expected to have an impact on the field 
evolution in isolated as well as binary pulsars. 

Multi-polar Structure - All of our and similar investigations have been 
based on an assumption of a pure dipolar model for the magnetic field. 
Though calculations for isolated neutron stars does not show any apprecia- 
ble change in multi-polar structures (Mitra, Konar & Bhattacharya 1999) 
- the situation would change in presence of accretion or a very strong mag- 
netic field (Geppert et al. 1999) due to the importance of the Hall term 
requiring further investigation. 

The Magnetar Question - Amongst some of the more recent develop- 
ments the magnetars pose a great challenge for the existing theories of field 
evolution since they require a very rapid field evolution in isolated neu- 
tron stars. Though some work have already been done in this area (Heyl &; 
Kulkarni 1998, Geppert et al. 1999) - more detailed investigation is needed. 
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